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Abstract 
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orders, and the quadrupole-quadrupole interaction arising at the 2.5PN level. In particular, the 
specific effect of non-linear memory is computed using a simplified model of binary evolution in the 
past. The "plus" and "cross" wave polarisations at the 2.5PN order are obtained in ready-to-use 
form, extending the 2PN results calculated earlier by Blanchet, Iyer, Will and Wiseman. 
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I. INTRODUCTION 

A. On gravitational-wave astronomy 



With the ongoing scientific runs of the Laser Interferometric Gravitational Wave Obser- 
vatory (LIGO) 1] and the recent commissioning of the French- Italian VIRGO detector j^, 
there is tremendous excitement among physicists and astronomers about the possibility of 
the first ever direct detection of gravitational waves. The Japanese detector TAMA300 is 
a. r eadv taking data Q and the British-Gennan joint venture GEO600 wi„ be operational 
soon |4j . These detectors will form a world wide network capable of coincident observations 
and joint data analysis. The design sensitivity of the three LIGO interferometers together 
with that of VIRGO and GEO600 is sufficient for the detection of GWs from powerful as- 
trophysical sources like binary black holes in the mass range of about 40 M@ □ □ □ . These 
detectors are ground based and will be sensitive to the high frequency regime of the spec- 
trum of GW, i.e. between 10 Hz - 10 kHz. The proposed space-borne detector LISA (Laser 
nterferometric Space Antenna) will probe the universe at low frequencies (1Hz - 10~ 4 Hz) 
||. Even if LIGO-1 will not detect any gravitational waves, it will be able to put upper 
bounds on the fluxes and populations of different types of sources. The recent discovery of a 
binary pulsar PSR J0737-3039 9] has brought double neutron star coalescence rate estimates 
within an astrophysically relevant regime for LIGO-1 and VIRGO. The estimated event rate 
for compact binary merger for LIGO-1 can now go as high as 1 event per 1.5yr Q]. LIGO-2 
with its ten fold improved sensitivity should definitely be able to detect gravitational waves 
from different astrophysical sources. 

GW signals can be broadly classified into four categories: inspiral, periodic, burst and 
stochastic. The strategies for data analysis vary for different types of signals. Determin- 
istic signals from inspiral and periodic sources will be probed by matched filtering (to be 
discussed later), whereas statistical signals such as stochastic sources will be searched for 
by cross-correlating pairs of interferometers, seeking correlations and statistical variations. 
Unmodelled signals like those from supernova explosion or gamma ray bursts (GRBs) can- 
not in principle be looked for by any parametric techniques. One needs to monitor power 
excesses in the frequency domain or have notable amplitude variations in time to detect 
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them. In all these studies, coincidence between multiple detectors improves the efficiency of 
detection significantly. 

The joint observation of astrophysical phenomena in electromagnetic as well as non- 
electromagnetic windows (such as GWs and neutrinos) will have a huge impact on obser- 
vational astronomy in the future. The current theories of GRB's, involve either the merger 
of a compact binary or a collapse. In either case, the GRBs are associated with the emis- 
sion of gravitational radiation [11]. Detection of a merger GW signal would require a lower 



signal-to-noise ratio, if it is coincident with a GRB 



121 ] . The data analysis techniques for 



such phenomena will include comparison of the correlated output of two gravitational wave 



detectors, one prior to the GRB and the other at times not associated with the GRB 
The spectral flux distribution of the neutrinos from GRBs associated with supernovae, can 
test for the possible time delays between the supernovae and the GRB down to much shorter 
time scales than can be resolved with only photons Q]. 

In the cosmo logical context, the primordial gravitational waves resulting from inflation, 
contribute to the polarisation of the cosmic microwave background radiation (CMBR). Their 
observations could provide us valuable information about the early universe as well as im- 
portant checks on the inflationary cosmological model la ll 71. 18| . The proposed Decihertz 
Interferometer Gravitational Wave Observatory (DECIGO) |l9|. could even probe the accel- 
eration of the universe, and hence provide an independent determination of the spatial curva- 
ture of the cosmological model. This is important as projects like SuperNova/Acceleration 
Probe (SNAP) [20] would also probe the acceleration of the universe by electromagnetic 
means. 

Pulsars can spin up by accreting matter from the neighbouring star. Recent observations 
point towards a possibility of GW emission balancing the accretion torque putting upper 



bounds on pulsar spin 



2l\ . If these gravitational waves are due to r-mode or CFS instability 



in the pulsar, as it is presently understood, they can be detected by the advanced LIGO 

23j by advanced 



detector 22]. The detection of GWs from instabilities in relativistic stars 
LIGO would be able to test models of stellar interiors. 

The "chirp" from the inspiral of two compact objects is one of the most plausible GW 
signals the ground based detectors would look for. The early inspiral will fall in the sensitivity 
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band of the space-based detectors, where as, the late inspiral will be a good candidate source 
for the ground-based detectors. Though these GWs are extremely weak and buried deep in 
the detector noise, the large number of precisely predictable cycles in the detector bandwidth 
would push the signal up to the level of detection. One can then use the technique of 
matched filtering first for the detection of GW and later (i.e. off-line) for the estimation of 
the parameters of the binary. In order to have a good detection, it is extremely important 
to cross-correlate the detector output with a number of copies of the theoretically predicted 
signal (corresponding to different signal parameters) which is as precise as possible, and 
which remains in accurate phase with the signal. This has made general relativistic modelling 



of the inspiral 



ine; compact 



binary 



^yyyyyy 



ICB) one of the most demanding requirements for GW 



data analysis jo, |2 



B. An overview of current calculations of radiation from ICBs 

Currently, post-Newtonian (PN) theory provides the most satisfactory description of 
the dynamics of ICBs and gravitational radiation emitted by them. Starting from the 
gravitational- wave generation formalism based on multipolar post-Minkowskian expansions 
(see the next Section), the gravitational waveform and energy flux at the 2PN order 1 were 
computed by Blanchet, Damour and Iyer jsj]]. This incorporated the tail contribution at 
1.5PN order both in the waveform and in the energy flux; the polarisation states correspond- 
ing to the 1.5PN waveform were calculated in Ref. 3| (note that some algebraic errors in 
this reference are corrected in |35[). The 2PN results have been independently obtained 



33, 



The associated 



using a direct integration of the relaxed Einstein field equation 
polarisation states (i.e. the "plus" and "cross" polarisation waveforms) were obtained in 
Ref. J3| • These works provided accurate theoretical templates which are currently used for 
data analysis in all the laser interferometric GW detectors like LIGO/VIRGO. Extending 
the wave-generation formalism, the 2.5PN term in the energy flux, which arises from a sub- 

n 

dominant tail effect, was added in Ref. |3i|. In the case of binaries moving in quasi-elliptical 
orbits, the instantaneous parts of the waveform, energy flux and angular momentum flux at 



1 As usual the nPN order refers either to the terms ~ l/c 2n in the equations of motion, with respect to the 
usual Newtonian acceleration, or in the radiation field, relatively to the standard quadrupolar waveform. 
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2PN order were computed by Gopakumar and Iyer j^. The polarisations of the waveform 
at this order (in the adiabatic approximation) has been obtained more recently 3sj], and the 
phasing of binaries in inspiralling eccentric orbits has also been discussed (39 1 . 

The extension of the gravitational wave generation formalism to third post-Newtonian 
order, and computation of the energy flux up to 3.5PN accuracy, was achieved in 4Q]. To this 
order, in addition to the "instantaneous" contributions, coming from relativistic corrections 
in the multipole moments of the source, the results include several effects of tails, and tails 
generated by tails. But, unlike at the 2PN or 2.5PN order, where the calculation is free from 
ambiguities, at the 3PN order the incompleteness of the Hadamard self-field regularisation 
leads to some undetermined constants in the mass quadrupole moment of point particle 
binaries (we comment more on this below). On the other hand, the computation of the 
binary's flux crucially requires the 3PN equations of motion (EOM). These were obtained 
earlier by two independent calculations, one based on the ADM Hamiltonian formalism of 



general relativity 



41 



42 



43, 



equations in harmonic coordinates 



44 1, the other on t 



re direct 3PN iteration of the Einstein field 
48| . Both approaches lead to an undetermined 



constant parameter at 3PN when using a Hadamard regularisation, but this constant has now 

nn nnn 

been fixed using a dimensional regularisation [49, 50]. An independent method |51l. 15 21 53], 



using surface integrals together with a strong field point particle limit, has yielded results for 
the 3PN EOM in agreement with those of the first two methods. In particular, the EOM so 
obtained are independent of any ambiguity parameter, and consistent with the end result of 
dimensional regularisation 4^], sjj. The conserved 3PN energy is thus uniquely determined, 
and then the 3.5PN energy flux, together with the usual energy balance argument leads to 
the expression to the evolution for orbital phase and frequency at 3.5PN order (54J]. 

Once the gravitational wave polarisations are available, one proceeds to use them to 
construct templates for gravitational-wave data analysis. Templates for detection generally 
use the so-called restricted waveform (RWF), where the phase is at the highest available 
PN order, but the amplitude is retained to be Newtonian, involving only the main signal 
harmonics at twice the orbital frequency. Though for detection (which will be done online) 



the RWF approximation may be enoug 



h B, 7, 25, 3, 



the complete waveform will be useful 



for the parameter estimation (off-line). Recently, studies using the complete waveform, 
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which includes the contribution from higher harmonics besides the dominant one, have 



57 



The complete 



shown that it may play a vital role in parameter estimation 
waveform carries information which can break the degeneracy of the model, and allow one 
to estimate the otherwise badly correlated parameters. In the case of a chirping neutron star 
binary, the masses of the individual stars can be separated because of the mass dependence 



of the higher harmonics 



56| . In the case of black hole binaries, whose frequencies are too 



low to be seen in the detector sensitivity window for long, higher harmonics compensate for 
he information lost when the signal does not last long enough to be apparent in the data 
5(| • An independent study [5^ about the angular resolution of the space-based LIS A-type 
gravitational wave detectors with a time domain 2PN waveform, showed the importance of 
including higher PN corrections to the wave amplitude in predicting the angular resolution 
of the elliptic-plane detector configuration. 

In the present work we provide the gravitational waveform from ICBs to still higher 
accuracy, namely 2.5PN, which should in consequence be useful for future improved studies 
in GW data analysis, for both LIGO-type and LISA-type detectors. We shall include in the 
2.5PN waveform instantaneous as well as "hereditary" terms, exactly as they are predicted 
by general relativity, completing therefore the 2.5PN generation problem for binaries moving 

n 

in quasi-circular orbits initiated in 36 . Using the waveform we next obtain the two "plus" 
and "cross" GW polarisations at 2.5PN extending the results of [35J. We shall verify that 
the 2.5PN wave form is in perfect agreement, in the test-mass limit for one of the bodies, 
with the result of linear black hole perturbations . 

The plan of this paper is as follows. In Section m we outline the post-Newtonian gener- 
ation formalism, based on multipolar post-Minkowskian (MPM) expansions and matching 
to a general post-Newtonian source. In Section we discuss the structure of the gravi- 
tational waveform up to the 2.5PN order for general sources. Then we present our list of 
results for the mass and current-type multipole moments required for the computation of the 
waveform in the case of point particle binaries. In Section ITVj we compute all the hereditary 
contributions relevant to the 2.5PN waveform, which are made at this order of tails and the 
non-linear memory integral. The results for the 2.5PN waveform are given in Section |3 
ready to be implemented as 2.5PN GW template in LIGO/VIRGO experiments. Section 
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13 contains also some short discussion, concerning the multipole moments and future issues 
related with the 3PN waveform. 



II. THE POST-NEWTONIAN WAVE GENERATION FORMALISM 

The wave generation formalism relates the gravitational waves observed at a detector 
in the far-zone of the source to the stress-energy tensor of the source. Successful wave- 
generation formalisms mix and match approximation techniques from currently available 
collections. These include post-Minkowskian (PM) methods, post-Newtonian (PN) methods, 
multipole (M) expansions and perturbations around curved backgrounds. A recent review 
6(j] discusses in detail the formalism we follow in the computation of the gravitational field; 
we summarise below the main features of this approach. This formalism has two independent 
aspects addressing two different problems. The first aspect, is the general method applicable 
to extended or fluid sources with compact support, based on the mixed PM and multipole 
expansion (we call it a MPM expansion), and matching to some PN source. The second 
aspect, is the application to point particle binaries modelling ICBs. 

A. The MPM expansion and matching to a post-Newtonian source 

To define the solution in the exterior of the source within the complete non-linear theory 



we follow Refs. 
Thorne 



61 



62 



63 



64 



65 



who built on earlier seminal work of Bonnor |67( and 
to set up the multipolar post-Minkowskian expansion. Starting from the general 
solution to the linearized Einstein's equations in the form of a multipolar expansion (valid 
in the external region), we perform a PM iteration and treat individually each multipolar 
piece at any PM order. In addition to terms evaluated at one retarded time, the expression 
for the gravitational field also contains terms integrated over the entire past "history" of 
the source or hereditary terms. For the external field, the general method is not limited 
a priori to PN sources. However, closed form expressions for the multipole moments can 
presently only be obtained for PN sources, because the exterior field may be connected 
to the inner field only if there exists an "overlapping" region where both the MPM and 
PN expansions are valid and can be matched together. For PN sources, this region always 
exists and is the exterior (r > a) near (r <C A) zone. After matching, it is found that 



the multipole moments have a non-compact support owing to the gravitational field stress- 
energy distributed everywhere up to spatial infinity To include correctly these contributions 
coming from infinity, the definition of the multipole moments involves a finite part operation, 
based on analytic continuation. This process is also equivalent to a Hadamard "partie finie" 
of the integrals at the bound at infinity. 



The formalism, notably the asymptotic matching procedure therein, 



am, nas Deen e> 

0, 2Q, 3, Q 



3een explored 



The 



in detail and extended in a systematic way to higher PN orders 
final result of this analysis is that, the physical post-Newtonian (slowly moving) source 
is characterized by six symmetric and trace free (STF) time-varying moments, denoted 
{II, J Li Wl, Xl, Yl, Zl}, 2 which are specified for each source in the form of functional of 
the formal PN expansion, up to any PN order, of the stress-energy pseudo-tensor of the 
material and gravitational fields [72]. These moments parametrize the linear approximation 
to the vacuum metric outside the source, which is the first approximation in the MPM 
algorithm. In the linearized gravity case t^ v reduces to the compact-support matter stress- 
energy tensor T^ v and the expressions match perfectly with those derived in Ref. 

Starting from the complete set of six STF source moments {II, Jl, Wl, X l , Y l , Z l }, 
for which general expressions can be given valid to any PN order, one can define a different 
set of only two "canonical 1 source moments, denoted {Ml, Sl}, such that the two sets 
of moments {II, • • • , Zl} and {Ml, Sl} are physically equivalent. By this we mean that 
they describe the same physical source, i.e. the two metrics, constructed respectively out 
of {II, ■ ■ ■ , Zl} and {Ml, Sl}, differ by a mere coordinate transformation (are isometric). 
However, the six general source moments {II, ■ ■ ■ , Zl} are closer rooted to the source 
because we know their expressions as integrals over t^ v . On the other hand, the canonical 
source moments {Ml, Sl} are also necessary because their use simplifies the calculation of 
the external non-linearities. In addition, their existence shows that any radiating isolated 

nn 

source is characterized by two and only two sets of time- varying multipole moments [61, 68J. 

The MPM formalism is valid all over the weak field region outside the source including 
the wave zone (up to future null infinity). It is defined in harmonic coordinates. The 



2 As usual L = ii%2 ■ ■ ■ %i denotes a multi-index made of I spatial indices (ranging from 1 to 3). The integer 
I is referred to as the multipolar order. 
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far zone expansion at Minkowskian future null infinity contains logarithms in the distance 
which are artefacts of the harmonic coordinates. One can define, step by step in the PM 
expansion, some radiative coordinates by a coordinate transformation so that the log-terms 
are eliminated and one recovers the standard (Bondi-type) radiative form of the metric, 
from which the radiative moments, denoted {Ul, Vl}, can be extracted in the usual way 68]. 
The wave generation formalism resulting from the exterior MPM field and matching to the 
PN source is able to take into account, in principle, any PN correction in both the source 
and radiative multipole moments. Nonlinearities in the external field are computed by a 
post-Minkowskian algorithm. This allows one to obtain the radiative multipole moments 
{Ul, Vl}, as some non-linear functional of the canonical moments {Ml, Sl}, and then of 
the actual source moments {II, ■ ■ ■ , Zl}. These relations between radiative and source 
moments include many non-linear multipole interactions as the source moments mix with 
each other as the waves propagate from the source to the detector. The dominant non-linear 
effect is due to the tails of wave, made of coupling between non-static moments and the total 
mass of the source, occurring at 1.5PN order (~ 1/c 3 ) relative to the leading quadrupole 



radiation 



631 ] . There is a corresponding tail effect in the equations of motion of the source, 



occuring at 1.5PN order relative to the leading 2.5PN radiation reaction, hence at 4PN 
order (~ 1/c 8 ) beyond the Newtonian acceleration [^J. At higher PN orders, there are 
different types of non-linear multipole interactions, that are responsible for the presence of 
some important hereditary (i.e. past-history dependent) contributions to the waveform and 
energy flux. 

A different wave-generation formalism from isolated sources, based on direct retarded 
integration of Einstein's equations in harmonic coordinates, is due to Will and Wiseman 



and provided major improvement and elucidation of earlier investigations in the same 
line 0, 1A | ■ This formalism is based on different source multipole moments (defined by 
integrals extending over the near zone only), together with a different scheme for computing 
the non-linearities in the external field. It has currently been completed up to the 2PN 
order. At the most general level, i.e. for any PN extended source and in principle at any 
PN order, the Will- Wiseman formalism is completely equivalent to the present formalism 
based on MPM expansions with asymptotic matching (see Section 5.3 in j6J| for the proof). 
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B. Application to compact binary systems 



This application represents the second aspect of our approach. To this end, in the first 
instance, the compact objects (neutron stars or black holes) are modelled as point particles 
represented by Dirac 5-functions. Indeed one can argue that for compact objects the effects 
of finite size and quadrupole distortion induced by tidal interactions are higher order in 
the PN approximation. However, the general formalism outlined in Section 111 Al is set 
up for a continuous (smooth) matter distribution, with continuous T^, and cannot be 
directly applied to point particles, since they lead to divergent integrals at the location of 
the particles, when T^ int particle is substituted into the source moments {II, ■ • • , Z L }. The 
calculation needs to be supplemented by a prescription for removing the infinite part of the 
integrals. Hadamard regularisation, based on Hadamard's notion of partie finie, is what 
we employ. This is our ansatz for applying a well-defined general "fluid" formalism to an 
initially ill-defined point-particle source. 

To summarise: A systematic analytical approximation scheme has been set up for the 
calculation of waveforms and associated quantities from point particles to the PN order 
required (or permitted by given resources). A technical cost is the need to handle 5-functions 
in a non-linear theory, which is dealt with the Hadamard regularisation scheme or a variant 
of it. However, we already mentioned that at the 3PN order, subtleties arise due notably 
;o the so called non-distributivity of the Hadamard partie finie, which resulted, as shown in 
4ol ]. in some "ambiguities" when computing the 3PN mass-type quadrupole moment, which 
could entirely be encoded into three undetermined numerical coefficients £, k, (. These 
combined into the unique quantity 9 = £ + 2n + £ in the 3PN energy flux for circular orbits. 

The latter ambiguities are the analogues of the undetermined parameters found in the 
binary's EOM at 3PN order, namely ujk and u s in the canonical ADM approach j^l, 42|, and 
A in the harmonic-coordinates formalism |45l l46j | . The parameter A is related to the "static" 
ambiguity cu s by A = — — ^fr, while the "kinetic" ambiguity cu k has been determined 
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to the value || (see 



—30; 

44 |47| for details). The presence of the static ambiguity uj s or, 



equivalently, A, is a consequence of the Hadamard regularisation scheme which happens to 
become physically incomplete at the 3PN order. Recently, Damour, Jaranowski and Schafer 
4g| proposed to use a better regularisation: dimensional regularisation. This led them to 
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a unique determination of u s , namely uo s = 0. More recently 



the application of di- 



mensional regularisation to the computation of the EOM in harmonic coordinates has led 
to the equivalent result for A, which is A = — The EOM are thus completely deter- 
mined to the 3PN order within both the ADM and harmonic-coordinates approaches using 
Hadamard regularisation supplemented by a crucial argument from dimensional regularisa- 
tion in order to fix the last parameter. 3 All the 3PN conserved quantities are determined in 
. A 3PN accurate center-of-mass has been constructed and used to reduce the 



Rcfs. 



44, 



conserved energy and angular momentum |4£ 

The numerical values of the radiation-field- related ambiguity coefficients £, k and £ in- 
troduced in (^J have also been determined 2j| using dimensional regularization, so that 
the PN corrections to phasing are completely determined to 3.5PN accuracy. However, as 
we shall work in the present paper on the 2.5PN waveform, i.e. one half order before 3PN, 
these ambiguities will not concern us here, and we shall find that no ambiguity shows up in 
any of our calculations based on Hadamard's regularisation. In fact, it can be shown that up 
to the 2.5PN order Hadamard's regularisation as we shall employ here gives the same result 
as would dimensional regularisation. The reason is that, in the source multipole moments 
up to this order, there are no logarithmic divergences occuring at the particles' locations, 
which correspond in d dimensions to poles in e = d — 3. 

In the present work we will require, for the computation of the time-derivatives of mul- 
tipole moments, the EOM for the case of circular orbits to 2.5PN accuracy. We denote the 
PN parameter in harmonic coordinates by 

7 = —. (2-1) 
rc 

where r = |x| is the binary's separation (x = yx — y2 is the vectorial separation between 
the particles), and m = mi + m 2 is the total mass of the binary system. We pose X\ — — , 
X 2 = —, and v = XiX 2 . Occasionally we also employ 5m = m 1 —m 2 so that — = X\ — X 2 . 

3 Note that both calculations are performed in the limit e — > 0, where the dimension of space is 

d = 3 + e. The principle is to add to the end results given by the Hadamard regularisation |4lLl42l l45l l46| 
the difference between the dimensional and Hadamard regularisations, which is specifically due to the 
poles ~ l/e and their associated finite part. 
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Then the 2.5PN binary's acceleration reads (v = Vi — V2) 4 



dv 
~dt 



2 32 G 3 m 3 u 
Vx- — v + 0(6) 



(2.2) 



5 c 5 r 4 

where the explicit 2.5PN term (~ 1/c 5 ) is the radiation reaction force in the harmonic 
coordinate system used here. The radiation reaction force plays an important role in our 
calculation of the waveform — it must be consistently included in all replacements of accel- 
erations at 2.5PN order (however the reaction force yields no contribution to the energy flux 



at 2.5PN order for circular orbits 



3Jj|). In Eq. ()2.2|) the orbital frequency uj = 2n/P (where 



P is the orbital period) is related to the binary's separation r in harmonic coordinates with 
2PN accuracy by [zf| 



Gin 



OJ 



-3 + u 



7 + 



41 „ 
6 + — v + v 2 
4 



7 2 + 0(6) 



(2.3) 



In the following we shall also need the inverse of Eq. ()2.3j) . i.e. 7 in terms of lu, which can 
conveniently be written into the form 



7 = x < 1 



x + 
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1 v 
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x 2 + 0(6) 



(2.4) 



in which we have introduced the gauge invariant frequency-dependent parameter 



' Gmuj\ 2 ^ 



(2.5) 



III. THE 2.5PN GRAVITATIONAL WAVEFORM 



A. Waveform as a functional of multipole moments 



In an appropriate radiative coordinate system = (cT,X l ), the transverse-traceless 
(TT) projection of the deviation of the metric of an isolated body from flat metric defines the 
asymptotic waveform hj^ (lower-case Latin indices take the values 1, 2, 3). The leading-order 
1/R part of hj^ (where R — |X| is the distance to the body) can be uniquely decomposed 



4 We systematically use the shorthand 0(n) to mean a small post- Newtonian remainder term of the order 
of 0{c- n ). 
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68| into its radiative multipole contributions introduced in Section |H] Furthermore, the 
PN order of the asymptotic waveform scales with the multipolar order I. Hence, at any PN 
order only a finite number of multipoles is required, and we have, with 2.5PN accuracy, 



TT 

km 



2G 
c^R 



'Pijkm 



+- 

c 
1 

+- 

& 

1 

+- 

& 

1 

+ ? 

1 

+ ^r 



1 4 

-N a Uij a + -£ a 6(iVj>A'& 



—NabUijab + -£afe(i^f)ac-^6c 



1 2 

gQ ^abcU-ljabc + ^j- £ ab(iVj} a cdNbcd 



^QQ-^abcdUijabcd + gg^a&fi *j 



6(i v j)acde-^bcde 



2520 



NabcdeUijabcde + 



210 



ab(iVj) ac def^bcdef 



+ 0(6) 



(3.1) 



The LVs and V^'s (with L = ij • • • a multi-index composed of / indices) appearing in 
the above waveform are respectively called the mass-type and the current-type radiative 
multipole moments (see discussion in Section 111 Aj) . They are functions of the retarded time 
Tr = T — R/c in radiative coordinates, Ul(Tr) and Vl(Tr). We denote by N = X/i? 
the unit vector pointing along the direction of the source located at distance R from the 
detector. A product of components of N = (Ni) i=1 ^ t 3 is generally denoted N L = NiNj 
The Levi-Civita antisymmetric symbol reads e a bi, such that £123 = +1. The operator Vijkm 
represents the usual TT algebraic projector which is given by 



'Pijkm 2 (Pik'Pjm + 'Pim'Pjk 7~*ij7~*km) > (3.2a) 

Vn = 6 id - . (3.2b) 

Given an orthonormal triad (N, P, Q), consisting of the radial direction N to the observer, 
and two unit polarisation vectors P and Q, transverse to the direction of propagation, we 
define the two "plus" and "cross" polarisation waveforms by 
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h + = ±(P i P j -Q i Q j )h% £ , (3.3a) 
h x = l -{P l Q 3 +Q l P J )h^ , (3.3b) 

in which the projector Vijkm present in front of the TT waveform may be omitted. 

In the case of circular binary systems we shall adopt for P the vector lying along the 

intersection of the orbital plane with the plane of the sky in the direction of the "ascending 

node" A/", i.e. the point at which the bodies cross the plane of the sky moving toward the 

detector, and Q = N x P. Following the convention of Ref. 35^ the unit vector joining 

the particle 2 to the particle 1, i.e. n = (yi — Y2)/r where r = |yi — y2 1 , is given by 

n = P cos + (Q cos i + N sin i) sin 0, where i denotes the orbit's inclination angle and <fi is 

the orbital phase, namely the angle between the ascending node and the direction of body 

one. 5 The unit direction of the velocity, i.e. X such that v = ruo\ (for circular orbits), is 

I I 

given by A = — P sin0 + (Qcosi + Nsinz) cos0. (See Fig. 7 and Eq. (7.4) of |33l|.) 

Using the MPM formalism, the radiative moments entering Eq. (|3.1|) can be expressed 
in terms of the source variables with sufficient accuracy, that is a fractional accuracy of 
(9(6) = 0(c~ 6 ) relative to the lowest-order quadrupolar waveform. For this approximation 
to be complete, one must compute: mass-type radiative quadrupole L/y with 2.5PN accu- 
racy; current-type radiative quadrupole Vy and mass-type radiative octupole with 2PN 
accuracy; mass-type hexadecapole Uiju and current-type octupole with 1.5PN precision; 
U ijMm and V ijki up to 1PN order; U ijk i mn , V ijMm at 0.5PN; and finally U ijk i mno , V ijMmn to 
Newtonian precision. The relations connecting the radiative moments Ul and Vl to the 
corresponding "canonical" moments Mr and Sl (see Section III Al for a short recall of their 
meaning) are given as follows 



For the mass-type moments we have 



Hiienis ivi T. c 
Q, 65, 3- 

U ij{ T R ) = M^(T B ) + ^ £ dV [in (Z*-I) + 11] Mjf(V) 

dVMfl{V)Mfl{V) 



G [ 2 rTR 



5| 7 / a<i\" J j> 

1 m (5) ^ - 5 ?^ (4) - 2 ^ {3) ™ {2) -L- 1 - i r< 

7 



+-MflM 3>a - -MflM^ a - -M^M^ a + -e ab<t M^ a S b 



5 The angle 6 in our convention differs by J from the same in Refs. |3f|,|38(. We follow here the convention 



of BIWW 35], that is related to the BDI one 0,0 by 0bdi = <?!>biww 
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Uijkm(T R ) 



+0(6) 



U m (T R ) = M%(T R ) + 



2GM 



Tr 



dV 



In 



+0(5), 



M%L(T R ) + | | 2 M £V 



r fl -v\ 97 

26/ + 60 



26 



30 



+1 [ R dVM^M^V) 



21 ^( 5i 
+0(4) , 



-— M { ^M km> - —MyLM, 



63 - K4) w(l) 
g - -<ij Ivl km> 



— M (3) M (2) 



(3.4a) 



(3.4b) 



(3.4c) 



where the brackets <> denote the symmetric-trace- free (STF) projection, while, for the 
necessary current-type moments, 

2GM 



Vij(T R ) = S\f{T R ) 
+0(5), 







7" 


J dV 




1 — oo 




+ 6. 



(3.5a) 



In 



T R -V 
2b 



5' 

+ 3 



ssliv) 



+±e ab<i M®M k>b - l -e ab<l M^M^ b - 2S <i M J (4) 
+0(4) . 



(3.5b) 



[The underlined index a means that it should be excluded from the STF projection.] For all 
the other needed moments we are allowed to simply write 



U L {T R ) 
V L (T R ) 



mI\t r )+o($, 

sf{T R )+0$). 



(3.6a) 
(3.6b) 



In the above formulas, M is the total ADM mass of the binary system, which agrees with 
the mass monopole moment. The M^s and Sls are the mass and current-type canonical 
source moments, and , denote their p-th time derivatives. 

The parameter b appearing in the logarithms of Eqs. (J3.4|) and (J3.5j) is a freely specifiable 
constant, having the dimension of time, entering the relation between the retarded time 
T R = T — R/c in radiative coordinates and the corresponding one t — p/c in harmonic 
coordinates (where p is the distance of the source in harmonic coordinates). More precisely 
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r ^ t _P_^£ ln (4). (3.7) 



we have 

T^ — i. — tL 

c c 3 \cb< 

The constant 6 can be chosen at will because it simply corresponds to a choice of the origin 
of radiative time with respect to harmonic time. 

As recalled in Section Hi Al the "canonical" moments Ml, Sl do not represent the best 
definition for what should be referred to as the "source" moments. This is why we now 
replace the M^'s and Sl's by adequate source multipole moments II, Jl, Wl, • • ■ , which 
admit closed-form expressions in terms of the source's stress-energy tensor. At the 2.5PN 
order it turns out that we need only to take into account a correction in the 2.5PN canonical 
maSS - tyP e quadrup0l e _ aad gl ven, , a center &am e, by (see Re fe . HQ, • 



AG 



0(7), (3.8) 
where iy is the source mass quadrupole, and where W denotes the "monopole" corresponding 



to the set of moments Wl- [We shall need W only at the Newtonian order where it will 
be given by (j3.18b|) ; see Section IIII CI for the expressions of all the source moments in the 
case of circular binary systems.] Note that a formula generalizing Eq. (J3.8j) to all PN orders 
(and multipole interactions) is not possible at present and needs to be investigated anew 
for specific cases. This is why it is more convenient to define the source moments to be II 
and Jl (and the other ones Wl, ■ ■ ■ , Z L as well, but in view of e.g. Eq. (|3.8|) these appear 
to be much less important than II, Jl) rather than Ml and Sl- For all the other moments 
needed here, besides the mass quadrupole (j3.8jl . we can write, with the required precision, 
that Ml agrees with the corresponding II and that similarly Sl agrees with Jl- Namely we 
always have 

M L = h + 0(5) , (3.9a) 
Sl = Jl + 0(5), (3.9b) 

and we can neglect in the 2.5PN waveform all the remainders in ()3.9|) except for the case of 
Mij where the required result is provided by Eq. (|3.8jl . Thus, from now on, the waveform will 



6 The equation (11.7a) in 40] contains a sign error, but with no consequence for any of the results of that 
reference. The correct sign is reproduced here. 
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be considered a function of the "main" source moments II, Jl, and also, of the "auxiliary" 
moment W appearing in Eq. ([3.8)1 . 



B. Instantaneous versus hereditary contributions 

From Eqs. (|3.4|) - (|3.5|) it is clear that the radiative moments contain two types of terms, 
those which depend on the source moments at a single instant, namely the retarded time 
Tr = T — R/c, referred to as instantaneous terms, and the other ones which are sensitive to 
the entire "past history" of the system, i.e. which depend on all previous times (V < Tr), 
and are referred to as the hereditary terms. 

In this work, we find it convenient to further subclassify the instantaneous terms in the 
radiative moments into three types based on their structure. The leading instantaneous 
contribution to the radiative moment Ul (resp. Vl) from the source moment II (resp. Jl), 
is of the form I® or jf 1 . We refer to these as instantaneous contributions from the source 
moments and denote them by the subscript "inst(s)". Starting at 2.5PN order, additional 
instantaneous terms arise, of the form I^ Ikm or nf 1 ^ki m the expressions relating radiative 
moments to source moments [see Eqs. ()3.4|) - (j3.5J) ]. We call such additional terms, the 
instantaneous terms in the radiative moment and denote them by the subscript "inst(r)". 
Thirdly, from Eq. (|3.8J) we see that the replacement of the canonical moments by the source 
moments induces also some new terms at 2.5PN level, of the form I^ W^, we shall call the 
instantaneous terms in the canonical moment denoted by the subscript "inst(c)". 

At the 1.5PN order, the hereditary terms are due to the interaction of the mass quadrupole 
moment with mass monopole (ADM mass M) and cause the effect of wave tails [63] . Physi- 
cally, this effect can be visualized as the backscattering of the linear waves (described by J^) 
off the constant spacetime curvature generated by the mass energy M. This can be viewed 
as a part of the gravitational field propagating inside the light cone (e.g. 64]). At higher 
PN orders there are similar tails due to the interaction between M and higher moments 1^, 
Jij, ■ ■ ■ . In addition, at the 3PN order (however negligible for the present study), there is 
an effect of tails generated by tails, because of the cubic interaction between the quadrupole 



moment and two mass monopoles, M x M x 4 



661 ] . The hereditary term arising at the 



2.5PN order in the radiative quadrupole ()3.4a)l is different in nature. It is made of the 
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quadrupole-quadrupole interaction, x I k i, and can physically be thought of as due to the 
re-radiation of the stress-energy tensor of the linearized quadrupolar gravitational waves. 



77 



78, 



It is responsible for the so-called "non-linear memory" or Christodoulou effect 
(investigated within the present approach in jf^j, 6sj]). So far, all these effects are taken into 
account in the calculation of the waveform up to 2PN order 3jj and in the energy flux up to 



3.5PN 



36 



40 



54] . The two different types of hereditary terms will be denoted by subscripts 
'tail" and "memory". 

Summarizing, with the present notation, the total 2.5PN waveform may be written as 



(7,TT\ , /. TT\ , ( h TT\ , (rTT\ , /. TT\ 

\ n km) inst ( s ) ^ V L km) inst ( T ) ^ V l km) hmtM ^ V W tail ^ V L km) 



t(r) 



(c) 



memory 



+ 0(6). 



(3.10) 



We give each of the above contributions explicitly. The instantaneous part of type (s) reads 
(h TT ) 

Vfcm/inst(s) 



2G 1 (2) 

ijkmi Iaa 



1 

+ - 

C 

1 

+ - 

1 

+-; 



+ - 



-Njff a + -e ab{i J^N b 



r(2) 



- at r( 4 ) 2c , T^°' AT 
Y2 ab ' ijab ~ l ~ 2 ^ j)ac bc 



r(3) 



' AT r(5) ~ - /> '< \- 



r(4) 



AT t( 6 ) AT 
ly abcd- l ija,bcd ' qf? tab ( it7 j')acde iVbccie 



(5) 



360 



36 



1 AT r(7) , J__ / t(6) n 
2520 bcde ijabcde "r 2]_g afe ( J J j)acdef^ bcde f 



where all the source moments are evaluated at the current time Tr. The type (r) is 



(3.11) 



TT 2G ^/ir(5) j _ 5.(4) _ 2 (; 

""Winstar) C 4 ^ ^ C 5 | 7 " ^ 1 a<i 1 j>a j^a 



(3) r (2) 



<i 1 j>a ^ t -ab<i 1 j>a u b 



+-N ab 



21 



~~ 7"( 5 ) r "-/''■' / 1 

^ *<ij*ab> g 1 <ij 1 ab> 



63 r(4) 

1 



102 r(3) r(2) 
5 1 <ij 1 ab> 



± F .r(B)r .rWrW _ 2/ rW 



. (3.12) 



Apart from two terms involving the source dipole moment Jj or angular momentum, these 
terms are made of quadrupole-quadrupole couplings coming from Uij, and in 
Eqs. (I3.4j) - (j3.5|) and computed in Ref. jf^ (though, using dimensional and parity argu- 
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ments, their structure can easily be written down, the computation of the numerical coeffi- 
cients in front of each inst(r) term needs a detailed study). The inst(c) terms refer to the 
instantaneous terms in the "canonical" moment and can be written down as 



{p'km) 



inst(c) 



2G G 

c4R r ijkmc5 < 



W (A Ifj + W (s)j(i) _ HfOO/W _ W (DjP)' 



, (3.13) 



y y tj 

where W is the particular "monopole" moment introduced in (|3.8jl . Both the inst(r) and 
inst(c) terms represent new features of the 2.5PN waveform. Concerning hereditary parts, 
we have the tail integrals (dominantly 1.5PN) which read as 



(Afem) 



tail 



2G 
c*R 



V, 



2GM 



ijkm 

N a 
+ — 

3c 

3c 

N ab 

12c 2 

Nbc 
2c 2 



In 
In 
In 
In 



Tr 

dV 

■oo 

Tr-V 

2b 
T R -V 

2b 
T R -V 

2b 
T R -V 
2b 



In 

97 
f 60 
7' 
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T R -V 

2b 

(4) 



11 

+ I2 



59 
f 30 



r (6) 

ijab 



(V) 



5 

+ 3 



and the non- linear memory integral (which is purely of order 2.5PN) given by 
■20 <* Tr 

' memory 



The latter expression is in complete agreement with the results of 



63, 



65 



(3.14) 



C. Source multipole moments required at the 2.5PN order 



Evidently the above formulas remain empty unless we feed them with the explicit expres- 
sions of the source multipole moments, essentially the mass- type II and current-type Jl, 
appropriate for a specific choice of matter model. In the present Section, we list the I^'s 
and Jl's needed for the 2.5PN accurate waveform in the case of point particles binaries in 
circular orbits. This is the extension of the list of moments given in Eqs. (4.4) of [3jJ for the 
computation of the 2PN accurate waveform. We do not give any details on this calculation 
because it follows exactly the same techniques as in Ref. |40j. 
Up to 2.5PN order in the waveform the mass moments are 
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Further, the current moments are given by 



J, 



z/m(X 2 -X 1 )STF ij \£ aM x : > a v 

1 



, , 13 4651 

+7 v 

1 9 252 



^mSTF^<! £fcafe x ai ^ b 



168 



1-3^ + 7 



,67 2 



+ 0(5). 
181 109 13 



90 



-z/ H ^ 

18 18 



7 r 2 

+ —^e kab x a v b ^(l-5u + 5u 2 
45 <r 



J, 



ijkl 



J ijklm 
Jijklmn 



V m 



(3.16a) 



(3.16b) 



(3.16c) 



(3.16d) 
(3.16e) 
(3.16f) 



(3.17a) 



(3.17b) 



!/ m (X 2 - XO STF iifei { £ ^'*V [1 - 2u 

+7 (ll - lT + Tl" 2 ) } + Ti ^ vkb i 1 - ^ + 3 ^ } 

+ 0(3) , (3.17c) 
STFy Wm {^x^V (l - 5i/ + 5z/ 2 ) } + 0(2) , (3. 17d) 



!/ m (X 2 - Xx)(l - 4z/ + 3z/ 2 ) STF yWmn {£ M ^ ffl3Hm !i ft } + 0(1) . (3.17e) 



[We recall that X x = ^, X 2 = , and v = XiX 2 ; the PN parameter 7 is defined by fl2~TJ) : 
the STF projection is mentioned explicitly in front of each term.] 

In addition, the current dipole J; in (|3.12j) is the constant binary's total angular mo- 



20 



mentum which needs to be given only at Newtonian order, and we need also to give the 
monopolar moment W which appears inside the inst(c) terms of (j3.13|) and comes from the 
relation (|3.8|) between canonical and source quadrupoles. We have 

Ji = ume iab x a v b + 0{2) , (3.18a) 
W = ^z/mx.v + 0(2). (3.18b) 

With all the latter source moments valid for a specific matter system (compact binary in 
circular orbit) the gravitational waveform is fully specified up to the 2.5PN order. 

IV. HEREDITARY TERMS AT THE 2.5PN ORDER 

We now come to the computation of the hereditary terms, i.e. made of integrals extending 
on all the past history of the non- stationary source, at any time V from — oo in the past up 
to Tr = T — R/c. In the following we shall refer to Tr as the current time — the one at 
which the observation of the radiation field occurs. As we have seen in Section IIIIB1 at the 
2.5PN order the waveform contains two types of hereditary terms: tail integrals, given by 
Eq. (|3.14j) . and the non- linear memory integral ()3.15j) . 7 

Evidently, in order to evaluate the hereditary terms, one must take into account the fact 
that the binary's orbit will have evolved, by gravitational radiation reaction, from early time 
on. However we shall show, following Refs. [3, &l] , that the tails can basically be computed 
using the binary's current dynamics at time Tr, i.e. a circular orbit travelled at the current 
orbital frequency oj(Tr) (this will be true modulo negligible 4PN terms in the waveform). 
Concerning the memory integral one definitely needs taking into account a model of binary 
evolution in the past. 

A. Model for adiabatic inspiral in the remote past 

In this paper we adopt a simplified model of binary's past evolution in which the orbit 
decays adiabatically because of gravitational radiation damping according to the standard 



7 Notice that in the gravitational- wave flux (in contrast to the waveform) , the non- linear memory integral 
is instantaneous — it is made of a simple time anti-derivative and the flux depends on the time-derivative 
of the waveform. 
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quadrupolar (i.e. Newtonian) approximation. We shall justify later that such a model is 
sufficient for our purpose — because we shall have to take into account the PN corrections 
in the tails only at the current epoch. The orbit will be assumed to remain circular, apart 
from the gradual inspiral, at any time V < Tr. We shall ignore any astrophysical (non- 
gravitational) processes such as the binary's formation by capture process in some dense 
stellar cluster, the successive supernova explosions and associated core collapses leading to 
the formation of the two compact objects, etc. 

Let us recall the expressions of the binary's orbital parameters as explicit functions of 
time V in the quadrupolar circular-orbit approximation |81]. The orbital separation r(V) 
evolves according to a power law, namely 

J m' J i> 



r(V) = 4 



G 3 r" 6: 



5 c 5 



■,1/4 



{Tc - V) 



(4.1) 



where T c denotes the coalescence instant, at which the two bodies merge together and the 
orbital frequency formally tends to infinity. The factor 1/c 5 therein represents the 2.5PN 
order of radiation reaction. We assume that our current detection of the binary takes place 
before the coalescence instant, Tr < T c , in a regime where the binary inspiral is adiabatic 
and the approximation valid. 8 

The orbital frequency uj in this model (uj = 2n/P, where P is the orbital period), is related 
at any time to the orbital separation (|4.1|) by Kepler's (Newtonian) law Gm = r 3 uj 2 . [Again 
we shall justify later our use of a Newtonian model for the early-time inspiral.] Hence, 



£5/3 5/3 



-3/8 



(4.2) 



Instead of uj it is convenient to make use of the traditional frequency-related post-Newtonian 
parameter x = (G muj /c 3 ) 2//3 already considered in (|2.5jl . It is given by 



<V) = - 



v c 



5Gm 



{Tc - V) 



-1/4 



(4.3) 



The orbital phase = J to dt = J x 3 ^ 2 dt, namely the angle between the binary's separa- 
tion and the ascending node Af, reads as 



<f>(V) = <Pc~ 



V c 
5Gm 



1 5/8 



(T c - V) 



(4.4) 



The formal PN order of the time interval left till coalescence is the inverse of the order of radiation-reaction, 
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where C is the value of the phase at the coalescence instant. 

The latter expressions are inserted into the various hereditary terms, and integrated from 
— oo in the past up to now. In order to better understand the structure of the integrals, it is 
advisable to re-express the above quantities (j4.1|) - (j4.4j) in terms of their values at the current 
time Tr. A simple way to achieve this is to introduce, following [8iJ|, the new time- related 
variable 

V = T^fr , (4-5) 

— J-R 

and to make use of the power-law dependence in time of Eqs. ()4.1j) - (j4.4|) . This leads imme- 
diately, for the orbital radius r(V) and similarly for uj(V) and x(V), to 



r(V)=r(T R )(l+y) 1 /\ 
where t(Tr) refers to the current value of the radius. For the orbital phase we get 



4>{V) = 0(T, 
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;i + yf* - 1 
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(4.7) 



5Gm 

The latter form is however not the one we are looking for. Instead we want to make explicit 
the fact that the phase difference between Tr and some early time V will become larger 
when the inspiral rate gets slower, i.e. when the relative change of the orbital frequency uj 
in one corresponding period P becomes smaller. 

To this end we introduce a dimensionless "adiabatic parameter" associated with the 
inspiral rate at the current time Tr. This is properly defined as the ratio between the 
current period and the time left till coalescence. We adopt the definition 9 

1 



Z{Tr) = 



(4- 



(T c — Tr) uj{T r ) ■ 

The adiabatic parameter £ is of the order 2.5PN. Written in terms of the PN variable x 
defined by (|4.3jl it reads 



256 v 



R 



x 5/2 (T R ). 



Now Eq. ()4.7|) can be expressed with the help of £(Tr) in the more interesting form 



<P(V) = ^{Tr) - 



;i+2/) 5/8 -i 



(4.9) 



(4.10) 



We have £ = |w/u; 2 so our definition agrees with the actual relative frequency change oc lj/uj 2 in one 



period. It is also equivalent to the one adopted in £ = Cbs/ 71 " 
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which makes it clear that the phase difference A0 = (j)(T R ) — <p{V), which is 2n times the 
number of orbital cycles between V and T R , tends to infinity when £(Tr) — > on any 
"remote-past" time interval for which y is bounded from below, for instance y > 1. What 
is important is that (j4.10j) depends on the current value of the adiabatic parameter, so we 
shall be able to compute the hereditary integrals in the relevant limit where £(Tr) — > 0, 
appropriate to the current adiabatic regime. Notice that, at recent time, when V — > T R or 
equivalently y — > 0, we have 

<f>(V) = <j>(T R ) - J— + {y 2 ) , (4.11) 
which is of course the same as the Taylor expansion 

4>{V) = <P(T R ) - (T R - V)u{T R ) + O [(T R - V) 2 ] . (4.12) 

B. The nonlinear memory integral 

We tackle the computation of the novel hereditary term at the 2.5PN order, namely the 
non- linear memory integral given by ()3.15|) . As we shall see the computation boils down 
to the evaluation of only two "elementary" hereditary integrals, below denoted I(T R ) and 
J{T R ). A third type of elementary integral, K(T R ), will be necessary to compute the tail 
integrals in Section II V CI 

The two wave polarisations corresponding to Eq. (j3.15|) . calculated with our conventions 
and notation explained after (|3.3|) . are readily obtained from the Newtonian approximation 
to the quadrupole moment [first term in (J3.16a)) ]. and cast into the form 

v /m . 2G G A m 5 v 2 sin 2 i [ Tr dV f 12 2 „ 
(h+) (T R ) = — / - - <^ 1 sin 2 i 

V + 7 memory V HJ ^ R g J ^ | 5 15 

+ ( TE ~ TE sin2 1 ) C0S ^ V ^ } ' ( 4 - 13a ) 

r, ^ x 2G G 4 m 5 u 2 sin 2 i [ Tr dV f4cosi . r „ 1/rr .,} , ^ , 

(^)_ y (T fi ) = m 3 ^ {— smmV))} ■ (4.13b) 

In our model of binary evolution the radius of the orbit, r(V), and orbital phase, <f)(V), are 
given by Eqs. (j4.1j) and ()4.4j) . at any time such that V <T R < T c . Because r(V) oc (T c — V) 1 ^ 
we see that the integrals in (J4.13|) are perfectly well-defined, and in fact absolutely convergent 
at the bound V — > — oo. There are two distinct types of terms in (|4.13|) . A term, present 
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only in the plus polarisation (|4.13a|) . is independent of the orbital phase 0, and given by a 
steadily varying function of time, having an amplitude increasing like some power law but 
without any oscillating behaviour. This "steadily increasing" term is specifically responsible 
for the memory effect. The other terms, present in both polarisations, oscillate with time 
like some sine or cosine of the phase, in addition of having a steadily increasing maximal 
amplitude. 

Consider first the steadily growing, non-oscillating term. Its computation simply relies, 
as clear from (j4.13a|) . on the single elementary integral 



where we find convenient to factorize out an appropriate coefficient in order to make it 
dimensionless. The calculation of ()4.14|) is easily done directly, but it is useful to perform 
our change of variable ()4.5|) . as an exercise to prepare the treatment of the (somewhat less 
easy) oscillating terms. Thus, we write in a first stage 



The factor in front is best expressed in terms of the dimensionless PN parameter x(Tr), and 
of course the remaining integral is trivially integrated. We get 



With this result we obtain the steadily-increasing or memory term in Eq. (j4.13a|) . How- 
ever, as its name indicates, this term keeps a "memory" of the past activity of the system. 
As a test of the numerical influence of the binary's past history on Eq. (j4.16|) . let us suppose 
that the binary was created ex nihilo at some finite initial time To in a circular orbit state. 
This premise is not very realistic — we should more realistically assume e.g. a binary cap- 
ture in some stellar cluster and/or consider an initially eccentric orbit — but it should give 
an estimate of how sensitive is the memory term on initial conditions. In this crude model 
we have to consider the integral Io(Tr) extending from To up to Tr. We find that the ratio 
of Io(Tr) and our earlier model I(T R ) is 




(4.14) 




(4.15) 




(4.16) 




(4.17) 
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Let us choose the observation time Tr such that the binary is visible by the VIRGO/LIGO 
detectors. At the entry of the detectors' frequency bandwidth, say tu se i S mic — 30 Hz, we 
obtain T c - T R ~ 10 3 s and r(T R ) ~ 700 km in the case of two neutron stars (m = 2.8 M Q ). 
For a binary system initially formed on an orbit of the size of the Sun, r(To) ~ 10 6 km 
(corresponding to T c — To ~ 10 8 yr), we find that the fractional difference between our two 
models Iq and / amounts to about 10~ 3 . For an initial orbit of the size of a white dwarf, 
r(To) ~ 10 4 km, the fractional difference is of the order of 10% — rather large indeed. So 
we conclude that indeed the memory term in (j4.13a|) . depends rather severely on detailed 
assumptions concerning the past evolution of the binary system. We shall have to keep this 
feature in mind when we present our final results for this term. 

Turn next our attention to the phase- dependent, oscillating terms in Eqs. (|4.13J) . Clearly 
these terms are obtained once we know the elementary integral 

(Gm) 4 r TR e 4 ^ y ) 
J(Tr) ^Ij V <_. (4 , 8) 

Inserting (j4.6J) and (j4.10|) into it we are led to the form [which exactly parallels (j4.16|) ] 

j < t *> = 2L x(T -> )e, ' t ' TR> r jttw-^ 1 ^'"' 11 ■ < 4i9) 

We shall compute this integral in the form of an approximation series, valid in the adiabatic 
limit £(Tr) — > 0. The easiest way to obtain successive approximations is to integrate by 
parts. We obtain 

-+oo 







dy 4| [(1+^5/8-!] = £_ f _ 7 f 00 dy ^[ (1+y) 6/«_ 1 ] 

;i + 2/) 5/4 4z\ 8j (1 + y)"/" 



(4.20) 



A further integration by parts shows that the integral in the curly brackets of ()4.20|) is itself 
of order £, so we have the result 10 

r^" ff[M,,H 4{ i+ °«)}- 



10 By successive integration by parts one generates the asymptotic series (divergent for any value of £) 
r<T^-W<-»>-'--.] ~ - 8 g(5„-3,(5„- 8 ,... (1 2,(r)(|)", 
where we have used the standard notation ~ for equalities valid in the sense of asymptotic series. 
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A standard way to understand it is to remark that, when y is different from zero, the phase 
of the integrand of (|4.21j) oscillates very rapidly when £ — > 0, so the integral is made of a 
sum of alternatively positive and negative terms and is essentially zero. Consequently, the 
value of the integral is essentially given by the contribution due to the bound at y = 0, 
which can be approximated by 

dye-~ = ^. (4.22) 

Because £(Tr) is of order 2.5PN the result 1)4.21)1 is sufficient for the control of the 2.5PN 
waveform, thus our elementary integral reads, with the required precision, 



J(T R ) = x 7 ' 2 {T R ) -^—{l + O (£)} . (4.23) 



We find that the "oscillating" integral J(Tr) is an order 2.5PN smaller than the "steadily 
growing" or "memory" integral I(Tr). This can be interpreted by saying that the cumulative 
(secular) effect of the integration over the whole binary inspiral in I(Tr) is comparable to 
the inverse of the order of radiation reaction forces — a quite natural result. By contrast the 
oscillations in J, due to the sequence of orbital cycles in the entire life of the binary system, 
compensate more or less each other yielding a net result which is 2.5PN smaller than for I. 
Furthermore, the argument leading to the evaluation of ()4.22j) shows that J, contrarily to 
J, is quite insensitive to the details of the binary's past evolution. 

Substituting Eqs. (J4.16)) and 1)4.23)1 into (j4.13j) we finally obtain the hereditary memory- 
type contributions to the polarisation waveforms as 

fi \ IGmvx . 2 .\ 17 + cos 2 % v 5/2 2 . I 

(h+) = — — sm 2 i< x b/2 (l + cos^) sin (40) > , 

(4.24a) 

(^Lemory = ^ ~ ^ * { ~ ^ C °^} ' 

where of course all quantities correspond to the current time Tp. The pha se- independent 



term is the non-linear memory or Christodoulou effect 



63 



65 



77 



pnat 

fi 



781 . [79j in the case of 

inspiralling compact binaries. Our calculation, leading to a factor oc sin 2 z(17 + cos 2 i), 

n 

agrees with the result of Wiseman and Will 78J. The non-linear memory term stricto- 



11 The difference of a factor of —2 of between their and our result is here probably due to a different 
convention for the polarisation waveforms. 
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sensu affects the plus polarisation but not the cross polarisation, for which it is exactly zero. 
As we have said, it represents a part of the waveform whose amplitude grows with time, but 
which is nearly constant over one orbital period. It is therefore essentially a zero-frequency 
(DC) effect, which has rather poor observational consequences in the case of the LIGO- 
VIRGO detectors, whose frequency bandwidth is always limited from below by (^seismic > 0. 
In addition, we know that detecting and analyzing the ICBs relies essentially on monitoring 
the phase evolution, which in turn is determined by the total gravitational-wave flux. But 
we have already noticed that the non-linear memory integral (j3.15|) is instantaneous in the 
flux (and in fact it does not contribute to the phase in the circular orbit case 3^|)- ^ seems 
thus that the net cumulative "memory" change in the waveform of ICBs is hardly detectable. 
On the other hand, the frequency- dependent terms found in Eq. (14. 24 J) form an integral part 
of the 2.5PN waveform. 

An important comment is in order. As we have seen the memory effect, i.e. the DC 
term in Eq. ()4.24a|) . is "Newtonian" because the cumulative integration over the binary's 
past just compensates the formal 2.5PN order of the hereditary integral. Thus we expect 
that some formal PN terms strictly higher that 2.5PN will actually contribute to the 2.5PN 
waveform via a similar cumulative integration. For instance at the 3.5PN level there will 
be a memory- type integral in the radiative quadrupole moment Uij, which is quadratic in 
the mass octupole moment (of the symbolic form J a f,<i x Ij> a b)- After integration over the 
past using our model of adiabatic inspiral, we expect that the "steadily-growing" part of the 
integral should yield a DC contribution to the waveform at the relative 1PN order. In the 
present paper we do not consider such higher-order post-Newtonian DC contributions, and 
leave their computation to future work. 

C. Gravitational- wave tails 

The tails up to 2.5PN order are given by Eq. (|3.14j) . Because of the logarithmic kernel 
they involve, the tails are more complicated than simple time anti-derivatives, and they 
constitute a crucial part of both the waveform and the energy flux, and in particular of the 
orbital phase with important observational consequences (for a review see 6o|). 



The computation of tails reduces to the computation of a new type of "elementary" 
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integral, differing from J(Tr) given in (j4.18|) by the presence of an extra logarithmic factor 
in the integrand, and given by 12 

K{T R ) = ^g- dV e —-\J^-4-\. (4.25) 



c 7 7-oo r^(V) \T C -T E 
For convenience we have inserted into the logarithm of (J4.25)) the constant time scale T c — Tr 
instead of the normalization 26 more appropriate for tails [or, for instance, 26e _11//12 , see 
(I3.14J) ]. but we can do this at the price of adding another term which will be proportional 
to J(Tr) already computed in Section llVBl 

With the help of the y- variable ()4.5|) we transform the latter integral into 

K(T R ) = JL- X (T R ) J«V e-rffa[<'«)"*-'] . (4.26) 

256 v J (1 + y) 5 / 4 

Because of the factor In y we are not able to directly integrate by parts as we did to compute 
J(Tr). Instead we must split the integral into some "recent-past" contribution, say y e] 0, 1 [, 
and the "remote-past" one, y el 1, +oo [. In the remote-past integral, whose lower bound at 
y = 1 allows for differentiating the factor In y, we perform the integrations by parts in a way 
similar to (|4.2U|) . This leads to 



+oo 



dy \ny 

■ e 



f [a + v)^-i] = j (427) 



k (1 + 2/) 5/4 

which is found to be of the order of O (^ 2 ) instead of O (£) in Eq. (|4.21|) . This is a consequence 
of the \wy which is zero at the bound y = 1, and thus kills the all-integrated term. Hence 
we deduce from (I4.27j) that the contribution from the remote past in the tail integrals is in 
fact quite small. The details concerning the remote-past activity of the source are negligible 
when computing the tails. More precisely, because £ = 0(5), we can check that terms 
such as ()4.27j) do not contribute to the waveform before the 4PN order, so the tail integrals 



12 



Actually in order to describe the tails we should consider the more general integrals 

(Gm)P- 1 fTn m4 ,(v) fT R -V 
K n , p {T R ) = / dV —— In 



2p - 3 y_oo rP(V) \T C -T R 



The dominant tails at the 1.5PN order correspond to p — 4 and n = 2, the tails at the 2PN order have 
p = 9/2 and n = 1,3, and those at the 2.5PN order have p = 5 and n = 2,4. Here we deal with the 
particular case p — 5, n — 4 because the calculation parallels the one of I(T R ) and J{T R ) in Section llV Bl 
The other cases are treated in a similar way. 
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can be legitimately approximat ed, with 2.5PN accuracy, by their "recent-past" history (in 



agreement with the findings of 32, 80]). 

Now, in the recent-past integral, i.e. y g]0, l[, we are allowed to replace the integrand 
by its equivalent when y — > 0, modulo terms of the same magnitude as (|4.27|) . This fact has 



been proved rigorously in the Appendix B of 
but simply state the end result, which reads 



. Here we shall not reproduce the proof 



jT'^^e-W^-l . f\ toj , e -¥ + 0ft . to{) . (4.28) 

As we see, the remainder is O (£ 2 ln£), instead of being merely O (£ 2 ) as in Eq. ([4.27)1 . The 
integral in the R.H.S. of ()4.28|) gives the main contribution to the tail integral (the only 
one to be taken into account up to very high 4PN order). It is easily computed by using a 
standard formula, 13 



1 _ 4^ y ^ 

dy \nye~~ = — 
4t 



2* U 



O (e 2 ) , (4.29) 



where C = 0.577 • • ■ is the Euler constant. 
Finally our needed result is 



K(T R ) = X V\T R ) — {- -h, (— ) - C + K Infl} , (4.30) 

where it is crucial that all the binary's parameters are evaluated at the current time Tr. It 
is interesting to compare (|4.30J) with our earlier result for J(Tr) given by (J4.23J) . in order 
to see the effect of adding a logarithmic-type kernel into an oscillating, phase- dependent 
integral. Eq. (|4.30|) . together with its trivial extension to K nm is used for the computation 
of all the tails in the waveform at the 2.5PN order (and it could in fact be used up to the 
order 3.5PN included). Actually we still have to justify this because during the derivation of 
(|4.3()jl we employed a Newtonian model for the binary's inspiral in the past, e.g. we assumed 
the Kepler law Gm = r 3 ( V) uj 2 { V) at any time V < Tr. In the case of the memory term 



13 



For any real number e we have 



e f dy lnye ley + i f + — e ley = - J sgn(e) - i (in |e| + C) , 
Jo Ji V 2 

where sgn(e) and |e| denote the sign and absolute value of e. 
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()3.15|) this is okay because it needs to be evaluated with Newtonian accuracy. But in the 
case of tails, Eq. ([3. 14)1 . the dominant effect it at the 1.5PN order, so we have to take into 
account a 1PN relative correction in order to control the 2.5PN waveform. Nevertheless, our 
model of Newtonian inspiral in the past is compatible with taking into account 1PN effects, 
for the basic reason that for tails the past behaviour of the source is negligible, so the 1PN 
effects have only to be included into the current values of the parameters, i.e. x(T R ), <fi{T R ) 
and £{T R ), in Eq. (|4~3H|) . 

To see more precisely how this works, suppose that we want to replace the "Newtonian" 
inspiral (|4.1j) by the more accurate 1PN law, 



r(V) 



G i m i v 
5c 5 



(Tc ~ V) 



1/4 



+ 



Gmr] 



(4.31) 



where r] = — — (see e.g. |6Cj). Substituting the variable y we obtain the 1PN 
equivalent of (|4.6|) as 

r(V) = r(T R )(l + y)Wi + r/[(l + y)" 1 / 4 - l]x(T R )\ , (4.32) 

where we parametrized the 1PN correction term by means of the variable x evaluated at 
current time T R (we consistently neglect higher PN terms). For the orbital phase we get 



<kv) = m 



5 i(T R ) 

743 11. 



5/8 



1 + (X(T R ))+T (1+2/) 



3/8 



x(T R ) \ , (4.33) 



foil + 24 u - ^ e mser f those expressions into our basic integral 



where £ = - — ~ 
()4.25|) . and split it into recent past and remote past contributions. Exactly like in ()4.27|) 
we can prove that the remote past of that integral, for which y e]l,+oo[, is negligible - 
of the order of 0{^ 2 ). Next, in the recent-past integral, y e]0, l[, we expand at the 1PN 
order [i.e. x(T R ) — > 0], to obtain some 1PN correction term, with respect to the previous 
calculation, of the form 



dy hiy 



(l + y) a -l 



-W [d+s/) 5/8 -i] 



(4.34) 



where a can take the values — |, | or |. Now the point is that this integral, like the remote- 
past one, is also of the order of 0(£ 2 ) or, rather, 0(£ 2 ln£). Indeed, the new factor {l+y) a — 1 
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in the integrand of ()4.34|) is crucial in that it adds (after taking the equivalent when y — > 0) 
an extra factor y, and we have thus to treat the following equivalent, 

/ dyylnye-^T =0(£ 2 ln£) , (4.35) 
Jo 

which is smaller by a factor £ = 0(5) than the integral ()4.29|) . as easily seen by integrating 
by parts. This means that the order of magnitude of the correction induced by our more 
sophisticated 1PN model for inspiral in the past is negligible. In conclusion, even at the 
relative PN order, one can use Eq. ()4.3())1 for computing the tails, but of course the current 
values of the binary's orbital parameters x(Tr), 0(Tr) and £(Tr) must consistently include 
their relevant PN corrections. This is what we do in the present paper, following the compu- 



tation in Refs. 
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of the higher-order tails up to relative 2PN order (i.e. 3.5PN beyond 
quadrupolar radiation) . 

Finally our results for the 2.5PN-accurate tail terms are as follows. It is convenient, 



following 



3, 



to introduce, in place of the "natural" constant time-scale b entering the tails 



and defined by (|3.7J) . a constant frequency-scale u given by 



ii-r 

- ■ (4-36) 



Like b, it can be chosen at will, for instance to be equal, as suggested in [35], to the seismic 
cut-off frequency of some interferometric detector, ujq = ^seismic- Then we find that all the 
dependence of the tails in the logarithm of the orbital frequency, i.e. the terms involving 
In uj and coming from the logarithm present in the R.H.S. of (|4.29jl . can be factorized out, 
up to the 2.5PN order, in the way 

where all the dependence upon \u(uj/ujq) is given as indicated [i.e. the (A;+, x )taii's are indepen- 
dent of ln(u;/cG>o)]- In the above expression the factor 1 — |x comes from the relation between 
the total ADM mass M and the binary's rest mass m = mi + m 2 , namely M = m [l — |x] . 
Because we are computing the tails with 1PN relative precision, this means that the factor 
of ln(a;/o;o), namely dh +)X /d(j) and therefore also /i +jX itself, is given at the relative 1PN 
order. The existence of this structure implies an elegant formulation of the 2.5PN waveform 
in terms of a new phase variable ip given by Eq. (j5.6J) below. The phase ip was already 
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introduced in |35(, and we have shown here that it is also valid, interestingly enough, for 
tails at the relative 1PN order. The "main" tails contributions are then given, up to 2.5PN 
order, by 

2 Gm vx 



( /C +)tail 



c 2 R 



{_2 7rX 3/2 (1+c 2 ) 



COS. 



+S— x 2 [(ll + 7^ + 10(5 + c?)ln2)si 
40 m 



sin i 



+x 



5/2 



-5tt(5 + c 2 ) cos - 27[7 - 10 ln(3/2)] (1 + c 2 ) sin 30 



-135vr(l + c 2 )cos3 



| (19 + 9c 2 - 2c\ + i/(-16 + 14c 2 + 6c*)) cos 20 



+ - (-9 + 14c 2 + 7c? + i/(27 - 42c 2 - 21c?)) sin 20 





16tt 



l-c*)(l-3z/) cos 40 



(*x) 



2 Gmvx 



tail 



f — (1 - c?)(l - 3z/) (21 - 201n2) sin- 

1 

{-4irx 3/2 Ci sin 20 



3 s,- & 5m 



x 2 [(3 + 10 In 2) cos + 5vr sin . 



20 m 

-9 [7 - 10 ln(3/2)] cos 30 - 45vr sin 30] 



+x 



5/2 



2tt 

— Ci (13 + 4 s 2 + i/(2 - 12 s 2 )) sin 20 
o 

— Ci(l - 3z/) (-6 + 11 s 2 ) cos 20 
5 

32tt 



Cj Sj (1 — 3z/) sin^ 



3 

+—Ci s]{l - 3i/) (-21 + 20 In 2) cos 40 
15 



(4.38a) 



(4.38b) 



Here Cj and Sj denote the cosine and sine of the inclination angle i, and 5m = mi — m 2 is 
the mass difference (so that 5m/m = X\ — X 2 ). Up to the 2PN order, we have agreement 
with the results of [35]. 

V. RESULTS FOR THE 2.5PN POLARISATION WAVEFORMS 



A. Instantaneous terms at the 2.5PN order 



In Section 1111 CI we have given the list of source multipole moments needed to control 
the waveform at the 2.5PN order. In this Section we proceed to calculate the instantaneous 
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terms [of types (s), (r) and (c)] in the 2.5PN waveform of circular compact binaries. The 
first step towards it is the computation of time derivatives of different moments II, Jl (and 
also W) using the binary's EOM up to 2.5PN order as given in Eqs. (|2.2|) - (j2.4j) . These time 
derivatives are then contracted with the unit direction N, for insertion into Eqs. (|3.11|) . 
(13.12)1 and ()3.13)1 . We can write down the inst(s) waveform schematically as, 

2Gvm 



inst(s) 



c 4 R 



(i) 



+(x 2 - Xl ) ^ + + (x 2 -x l ) > + plf ) } . 



(5.1) 



The instantaneous terms up to 2PN order were already reported in Eqs. (4.5) of Ref. (31 1. 
We have reproduced them in our present computation. At the 2.5PN order, most of the terms 
vanish when the two masses are equal (Xi = X 2 ) like at the previous "odd approximations" 
0.5PN and 1.5PN, but there is also an extra contribution, denoted by pf/ 2 ^ in Eq. ()5.1|) . This 



64 Gmj_ n (t v j) w \ l [ c \ l ar i S es directly from the 2.5PN 



consists of two parts: {p\J ) 

1 V v / reac 

radiation-reaction term in the EOM given by Eq. (|2~3j) : and (pi/ ') 



' quad 



_192 v GmrL n (i v j) 



7 



which comes from the 2.5PN contribution in the mass quadrupole fl3.1fiaj) . We find 

1408 Gm 7 2 



(5/2) 
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(5.2) 



The other contributions follow from a long but straightforward computation starting from 
the multipole moments listed earlier, and read as 
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where we recall that n = x/r and the parameter 7 is defined by (|2.1jl . 

In addition we must compute the instantaneous (r) and (c) parts of the waveform for the 
compact binaries in circular orbits. These parts are purely of order 2.5PN. The inst(c) part 
is computed starting from the expression for W in Eq. (j3.18b|) . but it turns out to be zero 
for circular orbits. We find 

z/7 2 f 84 Gm 



(h TT\ _ 2Gurn 



5 r 



■(N.n) (N.v) n 13 + 28(N.n)(N.v)v*- 
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inst(c) 
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(5.4a) 
(5.4b) 



B. The complete plus and cross polarisations 



We have computed all the five different contributions to the waveform contained in 
Eq. 1)3. 10J) . This together with the results of provides a complete 2.5PN accurate wave- 
form for the circular orbit case. In this Section, from the 2.5PN waveform, we present the 



35 



result for the two gravitational-wave polarisations, extending a similar analysis at the 2PN 
order in Ref. 3^|. 

The polarisations corresponding to the instantaneous terms are computed using Eqs. (J5.3|) 
and (|5.4|) . while those corresponding to the hereditary terms where obtained in Eqs. ()4.24j) 
and (|4.37jl - (|4.38jl . As in the earlier work, these polarisations are represented in terms of the 
gauge invariant parameter x = {Gmu / c 3 ) 2 / 3 , where lo represents the orbital frequency of 
the circular orbit, accurate up to 2.5PN order. This requires the relation between 7 and x, 
which has already been given in Eq. ([2.4)1 . The final form of the 2.5PN polarisations may 
now be written as, 

n (5 - 5) 

In particular we shall recover the 2PN results of |35J). However, for the comparison we have 



to employ the same phase variable as in 



which means introducing an auxiliary phase 



variable ip, shifted away from the actual orbital phase <fi we have used up to now, by Eq. 

is a priori adequate up to only the 2PN 



(5) of |3j|. Furthermore, the phase ifi given in j 
order, but we have proved it to be also correct at the higher 2.5PN order with the mass in 
front of the log-term being the ADM mass. Indeed, the motivation for the shift — > ip is 
to "remove" all the logarithms of the frequency [i.e. \nu or, rather, ln(a;/a;o)] in the two 
polarisation waveforms and to absorb them into the definition of the new phase angle ip. As 
a result, the two polarisation waveforms, when expressed in terms of ip instead of 0, become 
substantially simpler. From Eq. ()4.37j) we see that if we re-express the waveform by means 



ij) = (j)- 2x 3/2 1 - -x In ( — ) , (5.6) 





In 


f UJ \ 


1 - r. 











of the phase |3 



we are able to move all the lnw-terms into the phase angle. Notice that the possibility of 
this move is interesting because it shows that in fact the lnw-terms, which were originally 
computed as some modification of the wave amplitude at orders 1.5PN, 2PN and 2.5PN, 
now appear as a modulation of the phase of the wave at the relative orders 4PN, 4.5PN 
and 5PN. The reason is that the lowest-order phase evolution is at the inverse of the order 
of radiation-reaction, i.e. c 5 = 0{— 5), so as usual there is a difference of 2.5PN order 
between amplitude and phase. This shows therefore that the modification of the phase in 
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()5.6|) is presently negligible (it is of the same order of magnitude as unknown 4PN terms 
in the orbital phase evolution when it is given as a function of time). It could be ignored 
in practice, but it is probably better to keep it as it stands into the definition of templates 
of ICBs. The phase shift ()5.6|) corresponds to some spreading of the different frequency 
components of the wave, i.e. the "wave packets" composing it, along the line of sight from 
the source to the detector (see [80] for a discussion). 

With this above choice of phase variable, the same as in 3s| , all terms up to 2PN match 
with those listed in Eqs. (3) and (4) of 35j, though we recast them in a slightly different 
form for our convenience to present the 2.5PN terms. We find, 
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Notice a difference with the results of [35J], in that we have included the specific effect of 
non-linear memory into the polarization waveform at the Newtonian order, c.f. the term 
proportional to sf(17 + c?) in given by ()5.7aj) above. This is consistent with the order 
of magnitude of this effect, calculated in Section II V Bl However, beware of the fact that the 
memory effect is rather sensitive on the details of the whole time-evolution of the binary 
prior the current detection, so the zero-frequency (DC) term we have included into ()5.7a|) 
may change depending on the binary's past history (see our discussion in Section flVBj) . 
Nevertheless, we feel that it is a good point to include the "Newtonian" non-linear memory 
effect, exactly as it is given in Eq. (|5.7a|) . for the detection and analysis of ICBs. 14 

The purely 2.5PN contributions, in the plus and cross polarisations, constitute, together 
with the memory term in (|5.7ajl . the final result of this paper. They read as, 
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1 We already remarked that we have not computed the DC terms possibly present in the higher-order 
harmonics of the 2.5PN waveform. 
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Note that the latter cross polarisation contains a zero-frequency term [first term in Eq. 
(|5. IP)) ], which comes from the inst(r) contribution given by ()5.4j) . We employ the same 
notation as in 3^| , except that c$ and Sj denote respectively cosine and sine of the inclination 
angle % (which is defined as the angle between the vector N, along the line of sight from the 
binary to the detector, and the normal to the orbital plane, chosen to be right handed with 
respect to the sense of motion, so that < i < rr). In particular the mass difference reads 
5m = mi — m 2 . Like in [35(, our results are in terms of the phase variable ip defined by ()5.6|) 
in function of the actual orbital phase <fi (namely the angle oriented in the sense of motion 
between the ascending node M and direction of body one — i.e. = mod 2tt when the 
two bodies lie along P). We have verified that the plus and cross polarizations (J5.9J) - ()5.10J) 
reduce in the limit v — > to the result of black hole perturbation theory as given in the 



Appendix B of Tagoshi and Sasaki 



(the phase variable used in 



differs from ours by 



^ TS = jj, + tt/2 + 2x 3/2 [ln 2 - 17/12] and we have ^ TS = rr - i) 
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Equations ()5.9j) and 10)1 . together with ()5.7)) - (j5.8j) . provide the 2.5PN accurate tem- 
plate for the ICBs moving on quasi-circular orbits, extending the results of 3| by half a PN 
order. They are complete except for the possible inclusion of memory-type (zero-frequency 
or DC) contributions in higher PN amplitudes (1PN and 2PN). These wave polarisations 
together with the phasing formula of Ref. |54| . i.e. the crucial time variation of the phase: 
<f>(t), constitute the currently best available templates for the data analysis of ICB for ground 
based as well as space-borne GW interferometers. 



15 We spotted a misprint in the Appendix B of 59], namely the sign of the harmonic coefficient (be- 
having / = 7, m = 3, and corresponding to the cross polarisation) should be changed, so that one should 
129 „„ m\/i(;7 i u;v,/m/„5 



read (7^3 = + 1025024 o cos(0)(167 + ...) sm(9)(v 5 cos(3V>) - ...) 
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C. Comments on the 3PN waveform 



In Section Qll CI we have given the list of source multipole moments needed to control the 
waveform at the 2.5PN order. The computation of the 3PN waveform obviously requires 
more accurate versions of these moments as well as new moments, which all together would 



constitute the basis of the computation of the 3PN waveform. Thus, even though the 3PN 
mass quadrupole moment and 3PN accurate EOM Q, Q, Q, Q Q, Q, HQ, Q, Q 
are available, our present level of accuracy is not sufficient enough to compute the 3PN 
waveform. 16 

The source multipole moments at 3PN order, would yield the control of the inst(s) part 
of the waveform, as well as the tail terms, but we have to consider also other types of con- 
tributions, which are not all under control. The main reason for the present incompleteness 
at 3PN order is that the instantaneous terms of type (r) and (c), generalizing (|3.12|) - (j3.13|) 
to the 3PN order, are not computed. 

Recall that the inst(r) terms are those whose sum constitutes the instantaneous part of 
the relationships between the radiative moments Ul, Vl, and the "canonical" ones Ml, Sl, 
see (J3.4|) - (j3.5|) . Though one can guess the structure of these terms at the 3PN order, using 
dimensional and parity arguments, the numerical coefficients in front of each of them require 
detailed (and generally tedious) work. For instance, in the 3PN waveform we shall need the 
radiative mass-type octupole moment at the 2.5PN order, and therefore we have to know 
what is the remainder term 0(5) in Eq. (|3.4b|) which we don't (such a calculation would 
notably entail controlling the quadratic interactions between one mass and one current 
quadrupole, M^- x Ski, and between one mass quadrupole and one octupole, M^- x M^ m ). 
And similarly for the radiative current-type quadrupole V^- given by ()3.5a|) . We have also 
to compute the remainder terms 0(3) in the corresponding expressions of Uijkim and Vijki- 

Concerning the inst(c) terms, which are the instantaneous terms coming from the differ- 
ence between the canonical moments Ml, Sl and the general source ones II, Jl, • • • [c.f. 
Eqs. (I3.8j) - (j3.9j) ]. it does not seem to be obvious to guess even their structure. The crucial 



16 We are speaking here of the 3PN waveform. The computation of the 3PN flux is less demanding, because 
each multipolar order brings in a new factor c~ 2 = 0(2) instead of 0(1) in the case of the waveform, 



which explains why it is possible to control it up to the 3.5PN order in [4fJ. 
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new input we would need at 3PN order concerns the relation between the canonical mass 
octupole Mijk and current quadrupole Sij to the corresponding source moments 1^ and 



36|. 



at 2.5PN order, using for instance an analysis similar to the one in 

Finally, at 3PN order we would have to extend the present computation of hereditary 
terms. In the case of quadratic tails, like in 1)3. 14|) . the computation would probably be 
straightforward (indeed we have seen in Section II V CI that the complications due to the 
influence of the model of adiabatic inspiral in the past appear only at the 4PN order), but 
we have also to take into account the tail-of-tail cubic contribution in the mass-quadrupole 
moment at 3PN order, given in Eq. (4.13) of Ref. Q. In addition the analysis should be 
extended to the non-linear memory terms. The complete 3PN waveform and polarisations 
can be computed only after all the points listed above are addressed. 
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